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We develop a manifestly gauge-covariant expansion and projection using the eigen-mode of the 
QCD Dirac operator© = y^D^^ . Applying this method to the Wilson loop and the Polyakov loop, 
we perform a direct analysis of the correlation between confinement and chiral symmetry breaking 
in SU(3) lattice QCD calculation on 6"*^ at j3=5.6 at the quenched level. Notably, the Wilson 
loop is found to obey the area law, and the slope parameter corresponding to the string tension 
or the confinement force is almost unchanged, even after removing the low-lying Dirac modes, 
which are responsible to chiral symmetry breaking. We find also that the Polyakov loop remains 
to be almost zero even without the low-lying Dirac modes, which indicates the Z3 -unbroken 
confinement phase. These results indicate that one-to-one correspondence does not hold between 
confinement and chiral symmetry breaking in QCD. 
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1. Introduction: relation between confinement and chiral symmetry breaking 

Quantum chromodynamics (QCD) exhibits interesting nonperturbative phenomena such as 
color confinement and chiral symmetry breaking [|I|] in the low-energy region. In particular, it is an 
important issue to investigate the correlation between confinement and chiral symmetry breaking 
[Q, H, 0]. However, their relation is not yet clarified directly from QCD, although the strong cor- 
relation between them has been suggested by the simultaneous phase transitions of deconfinement 
and chiral restoration in lattice QCD both at finite temperature ^ and in a small- volume box 

The close relation between confinement and chiral symmetry breaking has been also suggested 
in terms of the monopole degrees of freedom [§, |3|], which topologically appears in QCD by taking 
the maximally Abelian gauge ||6|, ^, ||]. For example, by removing the monopoles, confinement and 
chiral symmetry breaking are simultaneously lost in lattice QCD [||], as schematically shown in 
Fig. 1. This indicates an important role of the monopole to both confinement and chiral symmetry 
breaking, and these two nonperturbative QCD phenomena seem to be related via the monopole. 
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Figure 1: An illustration of the relevant role of monopoles to nonperturbative QCD. In the maximally 
Abelian gauge, QCD becomes Abelian-like due to the large off-diagonal gluon mass of about IGeV [^], 
and there appears a global network of the monopole current ^]. By the Hodge decomposition, the QCD 
system can be divided into the monopole part and the photon part. The monopole part has confinement [^], 
chiral symmetry breaking [|]] and instantons [101, while the photon part does not have all of them. 



However, as a possibility, removing the monopoles may be "too fatal" for most nonperturbative 
properties. If this is the case, nonperturbative QCD phenomena are simultaneously lost by their cut. 

In fact, if only the relevant ingredient of chiral symmetry breaking is carefully removed, how 
will be confinement? To get the answer, we perform a direct investigation between confinement 
and chiral symmetry breaking, using the Dirac-mode expansion and projection JTl]]. 

2. Gauge-invariant formalism of Dirac-mode expansion and projection 

In this paper, we develop a manifestly gauge-covariant expansion/projection of QCD operators 
such as the Wilson loop and the Polyakov loop, using the eigen-mode of the QCD Dirac operator 
JD = y^D^, and investigate the relation between confinement and chiral symmetry breaking [pl|]. 
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2.1 Eigen-mode of Dirac operator in lattice QCD 

In lattice QCD with spacing a, the Dirac operator H) = y^D^ is expressed with U^{x) as 



1 ^ 

= ^lLyii [Uii{x)5,+f,^, - , (2.1) 



with U-^{x) = UJi{x — jx). Adopting hermitian y-matrices 7^ = 7^ , is anti-hermitian and satisfies 
= —Dx,y The normalized eigen-state \n) of the Dirac operator B is introduced as 

E)\n)=iXn\n) (2.2) 



with A„ € R. Because of {y5,E>} = 0, the state 75|?i) is also an eigen-state of E> with the eigenvalue 
— The Dirac eigenfunction V^„(jc) = {x\n) obeys lD\j/n{x) = iXnYn{x), and its explicit form of 
the eigenvalue equation in lattice QCD is 

1 

E [^M (■^) ¥n{x + fi)-U-^{x)Yn{x-fL)]= a,, Yn (x). (2.3) 

The Dirac eigenfunction \j/n{x) can be numerically obtained in lattice QCD, besides a phase factor. 
According to U^{x) V {x)U^{x)V^x + p.), the gauge transformation of Ynix) is found to be 

\l/„{x)^V{x)Yn{x), (2.4) 

which is the same as that of the quark field. To be strict, for the Dirac eigenfunction, there can 
appear an irrelevant «-dependent global phase factor as e"''" 1^1 , according to the arbitrariness of the 
definition of Yn{x). 

From the Banks-Casher relation [12], the quark condensate {qq), the order parameter of chiral 
symmetry breaking, is given by the zero-eigenvalue density p (0) of the Dirac operator E): 

{qq) = - lim lim 7rp(0), (2.5) 

m^OV^oo 

where the spectral density p(A) is given by p(A) = yT.„{S{X — with space-time volume V. 
Thus, the low-lying Dirac modes can be regarded as the essential modes responsible to spontaneous 
chiral-symmetry breaking in QCD. 

2.2 Operator formalism in lattice QCD 

The recent analysis of QCD with the Fourier expansion of the gluon field quantitatively reveals 
that quark confinement originates from low-momentum gluons below about IGeV in both Landau 
and Coulomb gauges [13|. This method seems powerful but accompanies some gauge dependence. 
To keep the gauge symmetry manifestly, we take the "operator formalism" in lattice QCD [11]. 

We define the link-variable operator by the matrix element of 

{x\U^\y) = U^{x)5x+f^^y. (2.6) 
The Wilson-loop operator W is defined as the product of along a rectangular loop, 

N 

W = llU^, = U^,U^,---U^,,. (2.7) 

k=i 
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For arbitrary loops, one finds Yjk=i fik = 0. We note tiiat the functional trace of the Wilson-loop 
operator W is proportional to the ordinary vacuum expectation value (W) of the Wilson loop: 

TrW = tiY,{x\W\x)=tiY,{x\U^,,U^,---U^,Jx) 

= tr 52 {M\U^t\x2){x2\U^^\x3){x3\U^^\x4.)---{xM\U^Jxi) 

Xt,X2,---^N 

2 N-1 

= tiY,{x\U^i\x + fii){x + fii\U^,_\x+Y,h)---{x+ £ h\U^Jx) 

X k=l k=l 

2 N-1 

= £tr{[/^,(x)[/^,(x + Ai)t/M3(^+LAi)---^Miv(^+ L A^)} = (W^)-Trl. (2.8) 

-v k=l k=i 

Here, "Tr" denotes the functional trace, and "tr" the trace over SU(3) color index. 

The Dirac-mode matrix element of the link- variable operator can be expressed with y/„ (x) : 

{m\U\n) =Y,{m\x){x\U^\x + jX){x + jX\n) =Y,¥l{^)U^i{x)Yri{x + jX). (2.9) 

X X 

Although the total number of the matrix element is very huge, the matrix element is calculable and 



gauge invariant, apart from an irrelevant phase factor. Using the gauge transformation ( |2.4| ), we 
find the gauge transformation of the matrix element as [ pl| ] 

{m\U^,\n) = £i/A,],(x)[/^(x)i/A„(x + /x) 

^ Y,wi{^)V\x)-V{x)U^{x)v\x + ^)-V{x + ^)\lfn{x + ^) 

X 

= Y^'^i{^)U^l{x)^ifn{x + fl) = {m\0^\n). (2.10) 

X 

To be strict, there appears an n-dependent global phase factor, corresponding to the arbitrariness of 
the phase in the basis \n). However, this phase factor cancels as e^"^"e^'^" = 1 between \n) and {n\, 
and does not appear for QCD physical quantities including the Wilson loop. 

2.3 Dirac-mode expansion and projection 

From the completeness of the Dirac-mode basis, ^„ \n){n\ = 1, arbitrary operator O can be 
expanded in terms of the Dirac-mode basis \n) as 

6 = Y^Y^\n){n\d\m){m\, (2.11) 



which is the theoretical basis of the Dirac-mode expansion Jl IQ. Note here that this procedure is 
just the insertion of unity, and is of course mathematically correct. 

Based on this relation, the Dirac-mode expansion and projection can be defined. We define the 
projection operator P which restricts the Dirac-mode space, 

P=Y^\n){n\, (2.12) 

neA 
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where A denotes arbitrary set of Dirac modes. In P, the arbitrary phase cancels between \n) and {n \ . 
One finds P = P and P' = P. The typical projections are IR-cut and UV-cut of the Dirac modes: 

^IR= L \n){nl PijY= £ \n){n\. (2.13) 

Using the projection operator P, we define the Dirac -mode projected link-variable operator, 

Uj^ =PU^P= £ £ \m){m\U^\n){n\. (2.14) 

meA neA 

During this projection, there appears some nonlocality in general, but it would not be important 
for the argument of large-distance properties such as confinement. From the Wilson-loop operator 
W = Y[k=i we define the Dirac-mode projected Wilson-loop operator W = Y[k=i ^lik' ^'^'^ 
rewrite its functional trace in terms of the Dirac basis as [ pl| ] 



N 

Tr = Tr n = Tr t/^, t/^, • • • f/^^ = Tr PU^,PU^,P ■ ■ ■ PU^.P 

k=i 



= tr {ni\U^i\n2){n2\U^2\>^3)---{nN\U^^\ni), (2.15) 

ni,n2,---,nfjeA 

which is manifestly gauge invariant. Here, the arbitrary phase factor cancels between 1?!^:) and (?ijt|. 
Its gauge invariance is also numerically checked in the lattice QCD Monte Carlo calculation. 
From Tr W^{R, T) on the RxT rectangular loop, we define Dirac-mode projected potential, 

V''{R) = - lim -ln{TrlV^(/?,r)}. (2.16) 



On a periodic lattice of V = L-^ x Nt, we define the Dirac-mode projected Polyakov loop []11|]: 
(LP) = -^Tr f]?7;= i-Tr {U^f' = ^tr £ {n,\U,\n2){n2\U4n,) ■■■ {n^^U^ln,), (2.11) 

i=i Ml,..,nAi,eA 

which is also manifestly gauge-invariant. 



3. Analysis of confinement in terms of Dirac modes in QCD 

In this paper, we mainly consider the removal of low-lying Dirac modes, i.e., the IR-cut case. 
Using the Dirac-mode expansion and projection method, we calculate the IR-Dirac-mode-cut Wil- 
son loop Tr W^{R, T), the IR-cut inter-quark potential V^{R), and the IR-Dirac-mode-cut Polyakov 
loop {Lp)i-R in a gauge-invariant manner |[IT|]. Here, we can directly investigate the relation between 
chiral symmetry breaking and confinement as the area-law behavior of the Wilson loop, since the 
low-lying Dirac modes are responsible to chiral symmetry breaking. 

As a technical difficulty, we have to deal with huge dimensional matrices and their products. 
Actually, the total matrix dimension of {m\U^\n) is (Dirac-mode number)^. On the lattice, the 
Dirac-mode number is x Nc-x 4, which can be reduced to be x A'^^, using the Kogut-Susskind 
technique [||]. Even for the projected operator, where the Dirac-mode space is restricted, the matrix 
is generally still huge. At present, we use a small-size lattice in the actual lattice QCD calculation. 



5 



Dime-mode expansion for eonfinement and chiral symmetry breaking 



Hideo Suganuma 



We use SU(3) lattice QCD with the standard plaquette action at j8 = 5.6 (i.e., a ~ 0.25fm) on 
6"* at the quenched level. The periodic boundary condition is imposed for the gauge field. We show 
in Fig.2(a) the spectral density p(A) of the QCD Dirac operator E>. The chiral property of ID leads 
to p(— A) = p(A). Figure 2(b) is the IR-cut Dirac spectral density Pir(A) = p(A)0(|A| — Ajr) 
with the IR-cutoff Air = 0.5a"' ~ 0.4GeV. By removing the low-lying Dirac modes, the chiral 
condensate is largely reduced as {qq) Air/ (qq) — 0.02 around the physical region of ntq ~ 5MeV. 
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Figure 2: (a) The Dirac spectral density p(A) in lattice QCD at j3=5.6 and 6^. The volume V is multiplied, 
(b) The IR-cut Dirac spectral density piR (A ) = p (A ) 6 ( | A | - Air) with the IR-cutoff Air = 0.5a" ' ~ 0.4GeV. 



Figure 3 shows the IR-Dirac-mode-cut Wilson loop {W''{R,T)) = TrW^{R,T), the IR-cut 
inter-quark potential V^{R), and the IR-Dirac-mode-cut Polyakov loop (Lp)iR, after the removal 
of the low-lying Dirac modes. These Dirac-mode projected quantities are obtained in lattice QCD 
with the IR-cut of Pir(A) = p(A)0(|A| - Air) with the IR-cutoff Air = 0.5a"' ~ 0.4GeV. 
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Figure 3: The attice QCD results after the removal of low-lying Dirac modes [jlip, which gives Pir(A) = 
p(A)0(|A| - Air) with the IR-cutoff Air = O.Sa"' ~ 0.4GeV. (a) The IR-cut Wilson loop Tr W^{R,T) 
(circle) after removing the IR Dirac modes, plotted against RxT. The slope parameter (7^ is almost the 
same as that of the original Wilson loop (square), (b) The IR-cut inter-quark potential (circle), which is 
almost unchanged from the original one (square), apart from an irrelevant constant, (c) The scatter plot of 
the IR-Dirac-mode-cut Polyakov loop (Lp)ir: its zero-value indicates Z3-unbroken confinement phase. 



Remarkably, even after removing the coupling to the low-lying Dirac modes, which are re- 
sponsible to chiral symmetry breaking, the IR-Dirac-mode-cut Wilson loop is found to obey the 
area law as (W^(/?,r)) oc e^^''^^ , and the slope parameter corresponding to the string tension 
or the confinement force is almost unchanged as ~ a. Accordingly, as shown in Fig.3(b), the 
IR-cut inter-quark potential V^{R) is almost unchanged from the original one, apart from an irrel- 
evant constant. Also from Fig.3(c), we find that the IR-Dirac-mode-cut Polyakov loop is almost 
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zero, i.e., (Lp)iR ~ 0, which indicates Z3-unbroken confinement phase. In fact, quark confinement 
is kept in the absence of the low-lying Dirac modes or the essence of chiral symmetry breaking 



|11[ ]. This result seems consistent with Gattringer's formula [|4p and Lang-Schrock's result [ ]14[ ] 



We also investigate the UV-cut of Dirac modes in lattice QCD, and find that the confining 



force is almost unchanged by the UV-cut [11], which seems consistent with the lattice result of 



Synatschke-Wipf-Langfeld JTsI]. Furthermore, we examine "intermediate-cut" of Dirac modes, and 



obtain almost the same confining force [|11[]. Then, we conjecture that the "seed" of confinement is 
distributed not only in low-lying Dirac modes but also in a wider region of the Dirac-mode space. 

Our lattice QCD results suggest some independence between chiral symmetry breaking and 
color confinement, which may lead to richer phase structure in QCD. For example, the phase tran- 
sition point can be different between deconfinement and chiral restoration in the presence of strong 
electro-magnetic fields, because of their nontrivial effect on chiral symmetry [|l^]. 
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